Weak Values Technique for Velocity Measurements 
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In a recent letter, Brunner and Simon propose an interferometric scheme using imaginary weak 
values with a frequency-domain analysis to outperform standard interferometry in longitudinal phase 
shifts [N. Brunner and C. Simon, Phys. Rev. Lett 105 (2010)]. Here we demonstrate an interfero- 
metric scheme combined with a time-domain analysis to measure longitudinal velocities. The tech- 
nique employs the near-destructive interference of non-Fourier limited pulses, one Doppler shifted 
due to a moving mirror, in a Michelson interferometer. We achieve a velocity measurement of 400 
fm/s and show our estimator to be efficient by reaching its Cramer-Rao bound. 



Introduction- In information theory, the Cramer-Rao 
bound (CRB) [U [5] is the fundamental limit in the min- 
imum uncertainty for parameter estimation. Measure- 
ments of phase |3H5] , beam deflection [7] , pulse arrival 
time [8j, Doppler shift [21 ^0\ and velocity pIUl4] are 
all fundamentally bounded by a CRB. If a measurement 
technique reaches the CRB, its estimator is said to be 
efficient. 

Spurred by fundamental studies of quantum phenom- 
ena |15| and by developments in precision measure- 
ments H [3 [3 [HI, the field of weak values [T7H2T] has 
become a powerful tool for parameter estimation [22H25] . 
The precision measurements inspired by weak values are 
not necessarily new or quantum. For example Zernike's 
phase contrast imaging [55], awarded the 1953 Nobel 
prize, can be classified as a weak value technique. An 
important aspect of the weak-values framework is that 
it provides a methodology for mitigating technical noise 
and amplifying an effect in one domain that is techno- 
logically difficult to observe in the conjugate domain. 

Recently, Simon and Brunner showed that a weak- 
values technique allows us to observe a large spectral 
shift induced by a small temporal shift. They also showed 
their technique outperforms standard interferometry lim- 
ited by technical noise _5]. Here we consider the oppo- 
site regime, where a small spectral shift causes a large 
temporal shift. The spectral shift in our experiment is 
a Doppler frequency shift produced by a moving mir- 
ror. Using established interferometry to measure veloc- 
ities arriving at the CRB is difficult but achievable as 
seen in Ref [2] . Our technique is comparable to standard 
interferometry, but allows us to reach its' CRB with a 
relatively simple method in a regime where 1// noise 
typically dominates. In this letter we show a weak- value 
optical technique to measure sub pm/s velocities. The 
protocal reaches the predicted CRB, averting technical 
noise and experimental imperfections. 

Theoretical description- The protocol, shown in 
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Fig. [T] uses a non- Fourier limited Gaussian pulse (i.e., 
CT ^ coherence length of the laser where r is the length 
of the pulse). The pulse, with initial intensity profile 
Iin{t) = -?^oexp (— i^/2T^), is sent through a Michelson 
interferometer with a slowly moving mirror in one arm. 
The interferometer is tuned slightly off destructive inter- 
ference by an amount 2(j), such that the output signal 
takes the form 
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where k = Svr/A, x{t) ~ vt and v is the velocity of the 
mirror. Assuming kvr ^ (f), making a small angle ap- 
proximation of (j) and re-exponentiating the output in- 
tensity, we obtain 
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Near-destructive interference reduces the peak intensity 
of the pulse by a factor sin^ (j), which is the probability 
for a single photon passing through the interferometer to 
reach the detector. Importantly, a time shift in the peak 
output intensity, 5t = 2kvT'^ /4>, has been induced with 
respect to the input. The velocity v can be obtained from 
measurements of the time shift St. 

We can rewrite the time shift in Eq. ([2]) in terms of 
the spectral shift St = 2kvT^ /cj) — 2TTfciT /(j) where the 
spectral shift, fd = 2v/X, of the pulse is proportional to 
velocity v. Instead of a direct spectral measurement, we 
obtain the velocity by measuring the induced time shift 
of the non- Fourier limited pulses. The time shift is ampli- 
fied in the measurement of v which is accompanied by a 
decrease in the measured intensity. These two results are 
well-known properties of the interferometric weak value 
amplification technique. In fact, a full weak value de- 
scription, using coherent Fourier limited pulses, can be 
formulated obtaining an identical result to Eq. ([2]). In our 
case, the use of non-Fourier limited pulses allows us to 
produce large time shifts regardless of the laser linewidth. 
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FIG. 1. An optical modulator generates a non-Fourier limited 
Gaussian shaped pulse. We couple the pulse to a fiber and 
launch it to a Michelson interferometer where one mirror is 
moving with constant speed v. The interference is controlled 
by inducing a phase offset 2(j} with the piezoelectric mirror. 
Photons exiting the interferometer are coupled into a fiber 
(not shown) and the arrival time of single photons are mea- 
sured with an avalanche photo diode (APD) and a photon 
counting module. 



We now consider the fundamental limitations of our ve- 
locity measurement set by the CRB. The CRB is equal 
to the inverse of the Fisher information, the amount 
of information a random variable (arrival time of pho- 
tons) provides about a parameter of interest (velocity). 
Assume that N photons are sent through the interfer- 
ometer. We want to determine the shift 6t from the 
set of N sin <j) independent measurements of photon ar- 
rival times. Such measurements follow the distribution 
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The CRB, F ^, is the minimum variance [571 [2H] of an 
unbiased estimation of St. The sensitivity in the deter- 
mination of St is therefore bounded by A {St) > t /(j) ^/~N . 
The error in the estimation of v is then bounded by 
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Note that this minimum uncertainty is independent of 
the actual value of v measured. This also determines 
the smallest resolvable velocity, when the signal-to-noise 
ratio is unity. The signal-to-noise ratio is 
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Experiment. -V^e use a grating feedback laser with 
A ~ 780 nm. An acoustic optical modulator creates 
Gaussian pulses of length t which we couple into a fiber. 
We launch them through the 50:50 beam splitter (BS) of 
the interferometer. The piezoelectric actuated mirror is 
driven by a triangle function with frequency /„i and peak- 
to-peak voltage Vpp. The pulse length is smaller than half 
the oscillating mirror period, so that a pulse experiences 
a single, constant velocity. An opposite constant veloc- 
ity is observed for each sequential pulse because the sign 
depends on whether the mirror moves toward the BS, or 
recedes away (see Fig. [l]) . The piezoelectric response a 
is calibrated by varying the voltage to change the dark 
port to a bright port. The piezo response was found to be 
a w 27 pm/mV for a low frequency-voltage product. The 
arm lengths (beam splitter-mirror distances) are approx- 
imately 1 mm (not including the BS size) to ensure long 
term phase stability. Photon arrival times are recorded 
with an avalanche photon diode (APD) and a photon 
counting module (PicoQuant PicoHarp 300). The detec- 
tor collects arrival times with 350 ps resolution. 

To calibrate the experiment we record the number of 
detected photons entering the interferometer, N . Then, 
the piezo-driven mirror is biased near destructive inter- 
ference and fed a triangle signal. We calculate the mean 
and error of the arrival time of the N^ detected photons 
for each set of pulses. The mean of the Gaussian de- 
termines the time shift St from which the velocity is ex- 
tracted, and the angle (p ~ yjN^/N is calculated. Lastly, 
to reach the CRB, we attenuate the peak of the pulses 
to about a million photons a second. 

Results.- We present velocity measurements w as a 
function of the pulse width r for different amplitudes on 
the moving mirror in Fig. [2] The lines are the theoret- 
ical predictions, v = 2fmVppa, where 2/,„ = l/6r. The 
mirror voltages are Vpp = {105, 52.5, 26.25, 10.5} mV 
and angle is = 0.31 ± 0.02 rad. The results agree well 
with the theoretical predictions. The smallest measure- 
ment of velocity in Fig. [Slis w = 60±11 pm/s. The angle 
(j) = 0.31 might seem large; however, comparing the exact 
form in Eq ([I]), | sin{(f> + kvt) / sin((/))|, to the approxima- 
tion, I exp(fcwt/(/))|, shows a discrepancy less than 1% for 
the experimental parameters. 

The uncertainties of the measurements in Fig. [2] are 
plotted separately in Fig. [3] and compared to the CRB 
Eq. (|4]). The error matches the CRB, thus the estimator 
is efficient and no other estimator can produce smaller 
uncertainties. This technique did not require noise fil- 
ters or frequency locking to reach the fundamental un- 
certainty in the mean arrival time of the photons. In 
addition, the fluctuations in the post selection angle 
are negligible. Therefore, our velocity measurement is 
fundamentally bounded by its CRB. 

It is important to note our CRB is scaled by the max- 
imum number of detected photons N . The collection- 
detection efficiency is about 20% due to the 50:50 BS 
(not shown in FigfTl) located before the APD used for 
alignment of the dark port, the efficiencies of the APD 
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FIG. 2. (color online). The Doppler shift, fd, is plotted as 
a function of r = {1.67, 4.17, 16.7, 417 and 833} ms. The 
phase offset angle is </> = 0.31 radians. The points are the ex- 
perimental results and the lines are the theoretical predictions 
for different voltages. Signal-to-noise ratios are 54, 27.4, 14.7, 
and 5.7 for Vpp ={105, 52.5, 26.25, 10.5} mV respectively. 



and the fiber coupling. Our calculations do not take the 
collection-detection efficiency into account. 

The results show precise and accurate detection of ve- 
locity measurements in the pm/s range. Results from 
Fig. [2] show smaller velocities can be measured with 
longer pulses. 

Now we seek to achieve the smallest velocities without 
the concern of reaching the CRB. Consider the temporal 
shift, advance or delay, of the pulse exiting the interfer- 
ometer. Since the peak of the pulse is sufficient to detect 
a the shift, we require a small region around the peak 
to determine the shift. This allows the use of effectively 
large values of r without requiring long term interfero- 
metric stability. Since the pulses are non-Fourier trans- 
form limited it is not necessary to use an entire Gaussian 
pulse. We truncate the Gaussian pulse to a width of r, 
that is 2/,„ = 1/t. In other words, the light intensity 
into the interferometer never drops below the 88% of the 
peak intensity and there is 12% peak to peak intensity 
variation following the peak of the Gaussian profile. 

We show velocities in the sub pm/s range using trun- 
cated pulses in Table [11 The mirror frequency was set to 
10 niHz, which corresponds to r = 50 s, and data was 
taken for voltages peak to peak, Vpp — (2, 1, 0.5} mV, 
for the piezo driving the mirror. Data was collected in 
intervals of 10 minutes (due to drift instability in inten- 
sity), and 13 sets of data were taken for each voltage. 
We did a Gaussian fit for each 10 minute interval. The 
time shift and its error were found as the mean and stan- 
dard deviation respectfully of the 13 time shifts obtained. 
The time shift was in the 10 s of millisecond and cor- 
responds to small Doppler shifts in the mircoHz range. 
This leads to the best technical noise limited measure- 
ment of (400 ± 400) fm/s. Nevertheless, both accuracy 
and precision are lost due to numerically fitting the trun- 
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FIG. 3. Experimental error in Fig. bl as a function of ryN. 
The solid line is the CRB as in Eq. (|4]) for A^ s; 54 x 10*^ 
photons. Note there are no error bars. 



cated distributions. Note that the measurements are all 
relative velocities because of the oscillating mirror. In 
one period there would be two pulses each with opposite 
but equal speeds. 

The results remain consistent with the full Gaussian 
picture theory, Eq. (l2|, but not with the CRB theory 
in Eq. Q. Calculating the mean arrival time of the 
photons is not a good estimator of the time shift be- 
cause we lack the full Gaussian pulse profile. Therefore 
we numerically fit the data to a unnormalized function 
Aexp [— (t — (5t)^/2r^] , the shift 5t is extracted and the 
velocity, u, is backed out. 



V^pp [mV] 


(f> [±0.002 rad] 


h [mHz] 


V [pm/s] 


2.0 


0.275 


3.6 ±1.2 


1.4 ±0.5 


1.0 


0.276 


1.6 ±1.1 


0.6 ±0.4 


0.5 


0.279 


1±1 


0.4 ±0.4 



TABLE I. Results of the cut Gaussian profile with r = 50 
s and A*' ~ 66 X 10^. The collection-detection efficiency is 
about 20%. The error is from the statistics of numerically 
fitting each run. Integration time was about two hours worth 
of data. 

Conclusion- In this letter, we show using non-Fourier 
limited pulses and standard interferometry inspired by 
weak values, sub pm/s velocities can be measured. Us- 
ing a Michelson interferometer tuned near a dark port 
we measure velocities as low as 400 ± 400 fm/s. We ac- 
complished sub pm/s velocity detection by bypassing the 
technical noise that flood intensity detectors to reach the 
CRB. The uncertainty of the phase measurement is neg- 
ligible when compared to the uncertainty of v for our 
parameter values so our uncertainty is the fundamental 
limit. Finally the error in our measurement of v matches 
the predicted CRB making this estimator efficient and 
the ultimate limit in uncertainty for velocity measure- 
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